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Abstract. We study the moments of the partial transpose of the reduced density 
matrix of two intervals for the free massless Dirac fermion. By means of a direct 
calculation based on coherent state path integral, we find an analytic form for these 
moments in terms of the Riemann theta function. We show that the moments of arbitrary 
order are equal to the same quantities for the compactified boson at the self-dual point. 
These equalities imply the non trivial result that also the negativity of the free fermion 
and the self-dual boson are equal. 
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1. Introduction 

The study of the entanglement content of extended quantum systems became in recent 
times a subject of extremely large theoretical interest (see e.g. the references in [1] as 
reviews). While the bipartite entanglement for an extended system in a pure state is a 
well understood subject and it can be quantified by the so-called entanglement entropy 
(i.e. the von Neumann entropy of the reduced density matrix of one of the two parts), for a 
mixed state the quantification of the entanglement is a more subtle issue. The observation 
that the entanglement in a bipartite mixed state is related to the presence of negative 
eigenvalues in the partial transpose of the density matrix [2] has led to the introduction 
of the negativity [3] which subsequently has been shown to be an entanglement monotone 
[4], i.e. a good entanglement measure from a quantum information point of view. 
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Although the negativity is a “computable measure of entanglement” [3], its direct 
and explicit computation in a many body system is very cumbersome. This difficulty can 
be (at least partially) overcome by a replica approach based on the computation of the 
even moments of the partially transpose density matrix [5]. This recent approach has 
been already applied to the study of one-dimensional conformal field theories (CFT) in 
the ground state [5, 6, 7, 8], in thermal state [9, 10], and in non-equilibrium protocols 
[10, 11, 12, 13], as well as to topological systems [14, 15, 16]. Focusing on ID CFTs in 
the ground state, the negativity is explicitly known only for the simple (but non trivial) 
geometry of two adjacent intervals embedded in a larger system. For the very important 
case of the entanglement between two disjoint intervals only the limit of close intervals 
is explicitly known. For arbitrary distances between the intervals, the main difficulty is 
to find the analytic continuation of the even-integer moments to n —> 1, although these 
moments are analytically known in a few cases [6, 7] (and indeed numerical interpolation 
techniques [17] have been exploited to have a numerical prediction for the negativity [18]). 
These analytical studies have been paralleled by several numerical works such as in Refs. 
[19, 20, 21, 22, 23, 24, 25], 

The goal of this manuscript is to investigate the negativity for the one-dimensional 
CFT of a free fermionic model. The main result is a close analytical form for the 
moments of the partial transpose of two disjoint intervals for the massless free Dirac 
fermion reported in Eq. (55). The explicit form of the moments is exactly the same as for 
the compactified boson at the self-dual point obtained in Ref. [6]. 

The manuscript is organised as follows. In Sec. 2 we build the partial transpose of 
the fermionic density matrix using coherent state path integral. In Sec. 3 we provide the 
analytical form for the moments of the partial transpose of two disjoint intervals and we 
analyse it. Finally in Sec. 4 we draw our conclusions and discuss some open problems. In 
a series of four appendices we report a number of technical details. 

2. Partial transpose of the reduced density matrix for the free fermion 

In this section we provide a path integral formula for the partial transpose of the density 
matrix for a free fermionic field theory, after a brief review of the result of Eisler and 
Zimboras [26] for the partial transpose of the reduced density matrix of two disjoint 
blocks on the lattice. 

2.1. Review of the lattice results 

We start from the the tight binding model with Hamiltonian 

1 L 

H = g [ c l c i+i + c I+i c * ’ 

i =1 

where periodic boundary conditions are assumed. We only consider the model at half 
filling kF = 7 t / 2 . Since (1) is quadratic in the fermionic operators, it can be diagonalized 


( 1 ) 
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in momentum space. The scaling limit of this model is the massless Dirac fermion which 
is a CFT with c—l. The local Hilbert space TLj of a single site is two dimensional and we 
can choose a basis made by the two vectors corresponding to whether the fermion occurs 
(|1)) or not (|0». In this basis the operators Cj and cj act as the creation and annihilation 
operators, Cj |0) = cj |1) =0, while cj |0) = |1) and c 3 |1) = |0). The tight-binding model 
can be mapped by a Jordan-Wigner transformation into the XX spin chain 


H xx = 


L 

n 

3 = 1 


a J a i+l + fyj+l 


( 2 ) 


Although these models are mapped one into the other, since the Jordan-Wigner 
transformation between them is not local, the entanglement (both entropy and negativity) 
of two disjoint blocks are not equal, as pointed out already in the literature [27, 28]. 

We always consider the entire system to be in the ground state T) with density 
matrix p = |T) (T|. It is useful to introduce the following Majorana fermions [29] 


a i = c 3 + C J’> a i = i( c i - c j)> ( 3 ) 

which satisfy the anticommutation relations {a", af} = 2 8 a p5 rs . The single site Majorana 
operators can be also written as 


4 = IX 


a y = -X 




(4) 


The operators X" on the single site satisfy the algebra of the Pauli matrices, but at 
different sites they anticommute and so they are not proper spin operators and should 
not be confused with the cr“ in (2). For each site, we also need to define the following 
unitary operator 

Uj® — e 1 ^ = cos (a/2) I + isin (a/2) Xj, (5) 

whose action on the Majorana operators (4) can be obtained from the following relation 

U™ X$ U { Q k) = [4,6 + (1 - 4,6 COS a)] X$ + (sin a)e kU Xj, (6) 


where £ k u is the totally antisymmetric tensor such that e xyz = 1. 

In this manuscript we are interested in a subsystem A = A\ U A 2 made by two disjoint 
blocks of lattice sites. Denoting by B = Bi U B 2 the complementary set of sites, which 
is also made by two disjoint blocks, the reduced density matrix is pa — Trg |\k) (J/|. This 
operator is Gaussian and it can be written as [29] 

Pa = ^efThJ2 Wl2 ° l0 ^ w l2 = (0\0\)i (7) 

where O k (with k G {1,2}) is a generic product of Majorana operators in A k , namely 

x y 

O k = njeAfc( a j)^ bl ( a j) My] w ith /irh G {0,1}. The sum in (7) is performed over all possible 
combinations of /i[^. 

Let us consider the operator 0 2 and introduce p 2 = 'f2j & A 2 (l- l [j] + /Ay) the total 
number of Majorana operators in 0 2 and p y 2 = P[j] the number of a^s in 0 2 . The 
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transpose of 0 2 (which obviously coincides with the partial transpose with respect to A 2 
in this case) is given by 


where 


t(M2) = 


(_1 )^(m)(_ 1 )^o 2j 

(8) 

0 (/i 2 mod 4) G {0,1}, 

1 (ji 2 mod 4) G {2, 3}. 

(9) 


The factor (—l) r ^ 2 ) in (8) originates from a rearrangement of the a x,y ’s operators after 
the transposition, while the factor (— 1)^2 comes from the fact that ( a y ) T = —a y and 
{af) T = af for the Majorana operators occurring in 0 2 . This extra factor can be removed 
by a unitary transformation. Another transposition can be naturally defined, namely 



(—1)^0 2 , 


( 10 ) 


where the unitary is now a product of terms like (5) over all the sites and it changes 
the sign of the a y, s leaving the a x, s untouched. This is the definition introduced in [26] 
and we will adopt this convention throughout this manuscript. Thus, let us drop the hat 
in (10) and denote it simply by Off. 

Given a block C of contiguous sites, an important ingredient in our analysis is the 
following string of Majorana operators 


P c II-K- 

3&C 


( 11 ) 


which satisfies Pff 1 = Pc- 

The partial transpose of pa with respect to A 2 can be written as the following sum 
of two Gaussian operators [26] 

Pa = ] ~y~ P A + ~Y~ p M~PaPa 2 = (e-^pA + e lf Pa 2 PaPa 2 ), (12) 

where the construction of pA has been reviewed in App. A and Pa 2 is the string of 
Majorana operators (11) along A 2 . 

The computation of Tr(pJ) n through (12) provides an expression containing 2 n terms 
given by all the combinations of p 0 = pA and pi = Pa 2 PaPa 2 i which can be written as 


Tr = 


i — 

e 4 




pi,P2,-..Pn=0,l 


2 n / 2 


Tr 


n Ppk 


L k =1 


(13) 


This formula can be further simplihed noticing that the various terms in the sum are 
invariant under the exchange pt —y 1 — p t . Using this and reorganising the terms in the 
sum, we can write 


Tr(rf?)" = 5SZT E 2 


n/2 


COS 


n—1 

2 ^ pi - n 


2=1 


Tr 


n—1 

p° n P pk 

■ k =1 


(14) 


where the vector p has n—1 components equal to 0 or 1 and therefore the sum contains 
2 n l terms. 
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In this subsection, we briefly review the features of the fermionic coherent states [30] 
which are needed to build the path integral of and pA- Here we focus on a single site 
(indeed, the site index will be dropped in this subsection) and in the next subsection the 
natural extension to many sites will be considered. 

The coherent states for fermions are defined through the Grassmann anticommuting 
variables. If Q\ and 6b are real Grassmann variables, we have that Of = 0 for i 6 {1,2} 
and Q\Q 2 = —0 2 0\. Since 6 2 = 0, a function f{6) of the real Grassman variable can be 
written as f{6) — fo + fid. Given two real Grassmann variables one can build a complex 
Grassmann variable ( as follows 

C = -^(^i + i0 2 ), C* = ;^= (0i - i # 2 ) ■ (15) 

The integration over a complex Grassmann variable acts as a derivation; indeed 


dC*dC = o, / dC*dC C = 0, / dC*dC c = 0, / dC’dC C C* = 1- (16) 


The coherent states are defined as follows 


1C) = |o) - C |i) , (Cl = (o| + C* (i| • (17) 

Since ( commutes with |0) and anticommutes with c, c 1 and 11), it is straightforward to 
check that c |C) = C 10 and (C| — (Cl C*- Notice that the coherent states do not provide 

an orthonormal basis. A completeness relation and a formula for the trace of an operator 


O read respectively 

I = /d<*d< e-« K) (Cl, TVO = y’dC*dCe-«(-C|0|C). (18) 

Given the above rules, the matrix elements of the identity and of the operators in 
(4) on the coherent states (17) can be computed, finding that 

(C\v) = 1 + Cv = {v*\~ C ), (19) 

(C\a x \v) = C + V = (v*\a x \C), (20) 

(C \a y \v) = -i(C* -v) = - (^*KIC*>, (21) 

((\ia x a y \g) = 1-(*V = «n m \M'a v \ ~ C*> , (22) 

where the second rewriting will be useful in the following subsection. Since i a x a y |C) = 
|—(), we can bring (20) and (21) in the same form of (19) and (22): 

«|V|,,} = <i*KI - O = -i (VI | - O , (23) 

«| 0 'W = (VK|-O= Hr\U ( 2 /2 a«ui%\-C). (24) 


(z) 

Notice that the insertion of Ufj 2 and its hermitian conjugate in (19) and (22) has no 
effect. 


Towards entanglement negativity for a ID free fermion 
2.3. Partial transpose of the reduced density matrix 


7 


The coherent state |£(x)) for a lattice is the tensor product of single site coherent states, 
with x runnig along the whole system or the corresponding subsystem. In the following we 
consider a lattice system but the final formulas can be extended to a continuous spatial 
dimension in a straightforward way by interpreting the discrete sums as integrals and 
integrations over a discrete set of variables as path integrals. 

The density matrix of the whole system in the ground state is p — |T) (T| and its 
matrix element between two generic coherent states reads 


p(C(x),v(x)) = e Cv (CO)I’I') (^\v(x)) , (25) 

where (*r] = C*( x ) 7 l( x )i with x labelling the whole system and eis the 

normalization factor (see (19)). To obtain the reduced density matrix in A, one first 
separates the degrees of freedom in A and the ones in B and then traces over the latter 
ones. Denoting by ICa^a)) and \(b(xb)) the coherent states on A and B respectively, we 
have that |C(x)) = |Ca(za)) ® \(b{x B ))- Adopting the notation \((x)) = |Ca(za), Cb(tb)), 
the matrix element of Pa is given by 

Pa(Ca(x a ),Va(x a )) = c~ C aVA J D X * b D Xb e~ x *B XB (Ca, ~Xb\^) {^\va, Xb) , (26) 


where D\*b Dxb = Y\. Xb d X*b( x b ) d Xb(xb) and the minus sign comes from the trace over 
B , according to (18). In the continuum limit, the braket (T|t 7 a, Xb) is the fermionic path 
integral on the upper half plane where the boundary conditions 7]a{xa) and Xb(xb) are 
imposed in A and B respectively, just above the real axis. Analogously, ((a, ~ Xsl'k) is 
the path integral on the lower half plane. The trace over B is performed by setting the 
fields along B equal (but with opposite sign) and summing over all the configurations. 
The resulting path integral is over the whole plane with two open slits along Ai and A 2 , 
where the boundary conditions t/a and (a are imposed respectively along the lower and 
the upper edge of A (left panel of Fig. 1). 

Let us consider the partial transpose p T f with respect to A 2 . In App. A the lattice 
results of [26] that we need in our analysis are briefly reviewed. Remembering that 
the partial transposition acts only on operators in A 2 , from (7) we can write its matrix 
elements as follows 

(C(x)Ipa \v(x)) = ^77^ (Ci(®i)|Oi|t7i(xi)) (Ci{x 2 )\Ol\r] 2 {x 2 )), (27) 

where Xj G A,-, with j G {1, 2}. 

Focussing on the term corresponding to A 2 in (27), from (10), (19)-(24) one finds 
<C 2 (^)| 0 2 t |^ 2 (x 2 )) = (— ir<«> (ri‘ 2 (xi )I ( f%0 2 ul% I - < 2 *(x 2 )) , (28) 


(z) 

where the unitary map U_^ 2 acts on all sites. When the number p 2 of Majorana operators 
in A 2 is even, from (9) we have that (—l) 1 "^ 2 ) = F 2 and therefore 
(Uxi)\(Ol)„ m \ m (x 2 )) = (-1F (r,;(x 2 )\ C/W /2 (0 2 ),v,„ ul% I - Q(x 2 )) 

= W 2 (x 2 )\ r/i»> f/W /2 (o 2 )„«„ ui% uj»> I - c 2 *(^)>, 


(29) 

(30) 
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Figure 1 . Path integral representation of p a and pa in the coherent state basis. 


where in (30) the factor (—l)^ 2 has been removed through a second unitary transformation 
which sends aj —> — a * leaving the aj’s unchanged (we recall that U^/ 2 exchanges the 
aj’s with the aj’s). The expression (30) suggests us to introduce the following unitary 
operator acting on A 2 

U = U® U% = n exp (- i \ = n exp 

jeA 2 ' 

whose net effect is to send a? —» —of and a 1 '{ —> —of, for j e A 2 . 

3 3 3 3 7 u 

In a similar way, we can treat the case of odd /r 2 , for which (—l) 1 "^ 2 ) _ ( see 

(9)). Again, from (10), (19)-(24) one gets 

<&(z2)l(OP od >(x 2 )} = - - Q(xP) (32) 

= - i (ig(i 2 )| V 2 (0 2 )odd | - Q(x 2 )) . (33) 

Introducing the operator 0 2 through its matrix elements as follows 

(C 2 (x 2 ) 1 0 2 \rj 2 (x 2 )) = (v 2 (x 2 )\ V 2 0 2 Vj | - C 2 (z 2 )>, (34) 


. TV ( -7T 4. _-7T 

(e 4 c] + e *c. 


. (31) 


the expression (27) can be written as follows 

(C(x)\Pa\v{x)) = 7fTJ^(^2 Wl2 <Ci(®i)|Oi|t7i(®i)) (C 2 (z 2 )1 0 2 \rj 2 (x 2 )) (35) 

' even 

-i^wis (Ci(®i)|Oi| 77 i(a;i)) (C2(^2)| 0 2 1772(^2)) Y 

odd ' 

where in the first (second) sum the terms have an even (odd) number of fermionic 
operators in A 2 (see App. A). In (35) the structure pj = p even — ip 0 dd (see (68)) can 
be recognised and this observation leads us to identify the matrix element of pa on the 
coherent states 


(C(x)\p A \v(x)) = (Ci(xi), 772 (^ 2)1 V 2 PaVJ |/ 7 i(xi), -£(x 2 )), (36) 

and analogously 

(C(x)\Pa 2 PaPa 2 \v{x)) = (Ci(®i), -^ 2 (^ 2 )! V 2 p A V% \pi(xi), ( 2 (x 2 )) , (37) 

where from (22) we can read that the action of Pa 2 is to change the sign of ( 2 . A 
graphical representation of this path integral representation for pA is given in the right 
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panel of Fig. 1. Hence, the final expression for the the partial transpose in the coherent 
state basis can be written exactly like the lattice counterpart i.e. 


Pa (Ci v) 


1 

71 


e ^Pa(CuC2;ViiV 2) + e 1 *p A (Cii-C2lVh-V2) , 


(38) 


where the notation is such that ry , Q £ Aj. This explicit form of the partial transpose in 
the coherent state basis is the final and main result of this section. 

In App. B we employ the formalism of this section to check the identity Tr(/rJ) 2 = 
Tr p\ (which holds for any quantum system [5, 6]) for the free fermion. 


3. Traces of the partial transpose for the free fermion 

In this section we consider the moments Tr(p^ 2 ) n for the free fermion. After summarising 
some needed CFT results for the moments of pa (Sec. 3.1), using the path integral 
approach of the previous section, we derive the analytic formula for Tr(p]J ) n given by 
(55) which is the main result of this manuscript. In Sec. 3.3 we show that the moments 
for the free fermion are equal to the ones for the compact boson at the self-dual radius. 
Finally, in Sec. 3.4 we give some numerical checks of our results. 


3.1. Some CFT results for Trp^ 

For the case of a single interval of length £ embedded in a CFT on the infinite line, the 
moments of the reduced density matrix can be written as [31, 32, 33] 

/ p \ —c(n—l/n)/6 

Trp^ = c n l-j , (39) 


where c is the central charge and a the inverse of an ultraviolet cutoff (e.g. the lattice 
spacing). The prefactors c n are non universal constants (that however satisfy universal 
relations [34]). The simple universal dependence on the central charge can be understood 
because Tr p n A is the partition function on a surface of genus zero that can be mapped to 
the complex plane [32], Eq. (39) can be interpreted as the two-point function of some twist 
operators acting at the endpoints of the interval u and v [32, 35], i.e. Tr p n A = {T n (u)T n (v)). 
The twist holds T n behave like primary operators with scaling dimension 

a ” = tK” 7)' (4o) 


The knowledge of the moments Trp^ give access to the full spectrum of the reduced 
density matrix [36]. While c n is not universal, its value for the tight-binding model at 
half-filling is known exactly and it is given by [37] 


= T 


exp < i n 


d^ log 


r (j + jg) 


[ tanh 


7 TZ 


tanh 


7 xnz 


)] 


(41) 
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For the case of two disjoint intervals A = A 1 U A 2 = [wi,rq] U [u 2 ,v 2 ], by global 
conformal invariance, in the thermodynamic limit, Tr p n A can be written as (dropping 
hereafter the dependence on the UV cutoff a) 


Tr p n A = cl 


\U i - U 2 \\Vi - V 2 \ 


\ui -Vi\\u 2 -V 2 \\ui “ v 2 11 U 2 V\ | 


2A n 


Fn(x), 


(42) 


where x is the four-point ratio (for real Uj and Vj, x is real) 


(ill - «l)(«2 - «2> 

(ill - ua)(»i - uj) 


6 ( 0 , 1 ). 


(43) 


The function F n (x) is a universal function (after being normalized such that .F n (0) = 1) 
that encodes all the information about the operator spectrum of the CFT while c n is 
the same non-universal constant appearing in (39). The function T n has been studied in 
several papers [38, 39, 40, 41, 42, 43, 44, 45, 46, 28, 47, 48, 49, 50] (see [51, 44, 52] for the 
holographic viewpoint and [53] for higher dimensional conformal field theories). In the 
case of two disjoint intervals, Tr p n A is the partition function on a surface of genus n — 1 
which cannot be mapped to the complex plane. This surface is usually called 7 Z n . 

One of the most important examples of exactly known F n {x) is the free boson 
compactihed on a circle of radius r circle . In this case, the function F n (x) (parametrized in 
terms of r/ = 2r^ ircle ) is [38] 


<d(0\r)T ) @(0|r/?/) 
[0(O|r)] 2 


(44) 


where r is an (n — 1) x (n — 1) matrix (called period matrix) with elements [38] 


' 1,3 


2 ^-4 

i — > sin(7rfc/n) 
n k =1 


2 Fi(k/n, 1 — k/n ; 1; 1 — x) 
2 Fi(k/n , 1 - k/n ; 1; x) 


cos[27r(/c/n)(i — j)\. 


(45) 


We remark that, since x G (0,1), the period matrix t(x) is purely imaginary. O is the 
Riemann theta function [54, 55] 

G(z\M) = J2 e m rnt ' M ' m+2iri ' z , (46) 

meZ "- 1 


which is a function of the (n — 1) dimensional complex vector z and of the (n — 1) x (n— 1) 
matrix M which must be symmetric and with positive imaginary part. 

For the critical Ising model, the scaling function F n (x) is also known [39] 


F n {x) 


1 

2 n ~ 1 0(O|r) 



e 

S 


(0|r) 


(47) 


where the period matrix r is the same as in Eq. (45). In this case 0 is the Riemann theta 
function with characteristic defined as [54, 55] 


Q[e](z\M) 


^ ^ g i7r(m+e) t -M-(m+£)-l-27ri (m+e) t -(z+5) 
m S Z n_1 


e = 



(48) 
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Figure 2. Cycles of type b for 1Z n (left) and lZ n (right). The solid line represents the 
part of the cycle belonging to the i-th sheet, while the dashed curve is the remaining 
part of the cycle, which lies in the (* + l)-th sheet. As for the cycles of type a, which are 
the same for 1Z n and 1Z n , we refer to Fig. 8 of [46]. 


where 2 and M are analogous to the ones in (46), and e, 8 are vector with entries 0 and 1/2. 
The sum in (e, 5) in (47) is intended over all the 2 n_1 vectors e and 8 with these entries. 
The parity of (48) as function of z is given by the parity of the characteristic, which is 
the parity of the integer number 4s ■ <5. There are 2 2 i" _ d characteristics: 2 n-2 (2 n ~ 1 + 1) 
are even and 2 n-2 (2 n ~ 1 — 1) are odd. In our following analysis only the trivial vector 
z = 0 occurs and therefore we will adopt the shortcut notation: @[e](M) = @[e](0|M) 
and @(M) = @(0|M) when the characteristic is vanishing. 

In the computation of the partition function on higher genus Riemann surfaces, one 
has to properly choose a canonical homology basis (i.e. a set of 2 (n — 1) closed oriented 
curves on the surface, the a and b cycles, which satisfy some specific intersection rules) 
and a set of n — 1 holomorphic differentials. By integrating such differentials along the b 
cycles one gets the period matrix of the Riemann surface. For a genus g Riemann surface, 
the period matrix is a g x g complex symmetric matrix with positive definite imaginary 
part [56, 57]. We refer the reader to [46] for a detailed analysis about the canonical 
homology basis for lZ n . In particular, the canonical homology basis {a r , b r ] 1 ^ r ^ n— 1} 
corresponding to (45) has been discussed in Sec. 4 of [46] and we will adopt it throughout 
this manuscript. In the left panel of Fig. 2 we show the j-th b cycle, which belongs to the 
j-tli sheet and to the (j + l)-th sheet. Instead, the construction of aj (which intersects bj 
only once) is more involved and therefore we refer the interested reader to Fig. 8 of [46]. 

The Riemann theta function with characteristic (48) occurs in the computation of 
fermionic models on higher genus Riemann surfaces [56, 57]. The characteristic e specifies 
the set of boundary conditions along the a and b cycles of the canonical homology basis 
and this provides the so called spin structures of the model. The vector e is determined 
by the boundary conditions along the a cycles (e*, = 0 for antiperiodic b.c. around a*, and 
£k = 1/2 for periodic b.c.), while 8 is provided by the boundary conditions along the b 
cycles (8k = 0 for antiperiodic b.c. around bk and 8k = 1/2 for periodic b.c.). 

3.2. Moments of the partial transpose for the free fermionic field theory 

We are finally ready to derive the moments of the partial transpose Tr(pj) n . The path 
integral for f/f is given by (38) which is a sum of two different operators. The moments 












Towards entanglement negativity for a ID free fermion 


12 


are then given by the sum of 2 n terms that come from the expansion of the binomial. 
Actually, since there is a double degeneration of these terms, the sum is only on 2 n ~ 1 terms. 
Introducing, in analogy with the lattice computation, the notation po(Ch) = Pa((,v) and 
Pi(C, v) — Pa (Ci ) — ( 2 ; Vi> —V 2 ) = Pa 2 PaPa 2 i the 2 n_1 terms in the sum for the moment of 
order n can be written as 

/ n n —1 

Y[Dx* k DxkPo(-Xn,Xi) Y[p Pk (Xk,Xk+i), (49) 

k =1 k=l 

with pi = 0,1. Each of these 2 n_1 terms is a partition function of a free fermion on a 
Riemann surface of genus n — 1 in which antiperiodic or periodic boundary condition are 
imposed along the basis cycles. 

At this point, before deriving the final result, we should discuss the Riemann surface 
7 Z n on which these partition functions are defined. The surface lZ n is defined by the 
density matrix pA and has genus n — 1 [5,6], but it is different from the one defining Tip n A 
(denoted in the previous section as 7 Z n ). Only for n = 2 they are the same torus (their 
moduli are related by a modular transformation), but for n > 2 they are different. The 
properties of this Riemann surface are discussed in details in App. C. The period matrix 
t(x) of 7 Z n for x G (0,1) is given by [6] 

f(x) = t(x/(x - 1)) = 1Z + il, 1Z=^Q, (50) 

where the elements of r have been defined in (45) and the real and imaginary parts of r(x) 
are 1Z and X respectively. Here it is important to observe that Q is a very simple symmetric 
integer matrix: it has 2 along the principal diagonal, —1 along the secondary diagonals and 
0 for the remaining elements. In App. C.l we report the detailed derivation of this result. 
As for the cycles of lZ n providing the canonical homology basis {a r ,b r ; 1 ^ r ^ n — 1} 
which gives the period matrix (50), we find that a r is the same as a r (we remind that lZ n 
and lZ n differ only for the way to join the sheets along A 2 ), while the generic cycle b r is 
obtained by deforming the cycle b r as shown in Fig. 2. 

An important ingredient at this point is the operator Pa 2 ■ For an arbitrary interval 
C we can write 

P c = = (-l) Fc . (51) 

where Fc is the fermionic number operator in the interval C which was already introduced 
long ago [58]. This operator is located along the interval C and it changes the fermionic 
boundary conditions (from antiperiodic to periodic or viceversa) on a cycle whenever it 
crosses the curve C. In (38) for pyy we have that Pa 2 occurs both before and after pA- 
This corresponds to the insertion of the operators (— 1)^2 above and below the cut along 

A-2- ~ 

Each term (49) is a partition function on 7Z n with some specific boundary conditions 
along the a and b cycles and it can be expressed in terms of Riemann theta functions. 
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Explicitly, we have 

0(f(a;)) 


Tr 


PO Ppk 


= C„ 


X 


U*2 



(52) 


where 0 is the vector made by n — 1 zeros. In this formula we have still to fix the vector 
<5 in terms of p , which is done as follows. Eq. (49) is evaluated on the n-sheeted Riemann 
surface lZ n where the i-th sheet is associated to the p Pi . On the sheets associated to pi, 
two operators (— 1) Fa 2 must be placed above and below A 2 . Then, the spin structure 
e can be read off by counting how many times the cycles of the basis cross the curves 
A 2 . Since the cycles a* do not intersect A 2 at all, we have that e = 0, i.e. the boundary 
conditions for the fermion along all the cycles a r are antiperiodic. Instead, for b r this 
analysis is non trivial because it intersects A 2 on the r-tli sheet and on the (r + l)-th 
sheet, as one can see from the right panel of Fig. 2. If b r crosses these curves an even 
number of times, then 6 r = 0, otherwise 6 r = 1/2. It is not difficult to conclude that 


2Si = (pi + p i+ 1 ) mod 2, 


(53) 


whose inverse reads 


Pi = 


n —1 


j=i 


26j i mod 2 


l-(-l) 


2 YTfZl h 


(54) 


The simplest example of (52) is the term Trp^ (namely p — 0). This spin structure 
has antiperiodic boundary conditions along all the cycles, i.e. e = S = 0. For this term 
D n [ef = 1. 

Thus, Tr (pj)” can be written as a sum over all the allowed spin structures: 


Tr {p T lT=cl 

The coefficient r n (S ) is 


l-x\ 2A " 1 






0[e](f) 


0(f) 


e = 


r n (6) = 2 n//2 cos 


n—1 


i + E *- 1 ) 22 


j=i h 


1=1 


(55) 


(56) 


It can be seen that r n (S) e {—2 n / 2 , 0, 2 n / 2 } for even n and r n (S ) G {—2^ n_1 ^ 2 , 2® -1 ^ 2 } 
for odd n. 

The analytic expression given by (55) and (56) is the main result of this manuscript. 
When the size of the intervals is very small with respect to their distance (di,d 2 -C d. i.e. 
x 1), it is possible to expand (55) in powers of x, as shown in App. C.2 where we find 
the first non trivial term of this expansion. 

There is also a very interesting by-product of our analysis which is given by (52) 
providing a very deep technical insight. Indeed Eq. (52) shows also that each of the 2 n ~ 1 
terms in the sum over p in (If) has a well defined continuum limit which is the partition 
function of the free fermion on lZ n with a particular assignment of fermionic boundary 
conditions, i.e. always antiperiodic along all the a cycles, while the b.c. along the b cycles 
are specified by <5 (we recall, antiperiodic for 6i = 0 and periodic otherwise). 
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3.2.1. Dihedral symmetry. The Riemann surfaces lZ n and 7 Z n enjoy a dihedral symmetry 
Z„ x Z 2 , as already noticed in [45]. The symmetry Z n comes from the invariance under 
cyclic permutation of the n sheets and the symmetry Z 2 corresponds to take the sheets 
in the reversed order and to reflect all of them with respect to the real axis. The 
former symmetry comes from the fact that 7 Z n and lZ n are obtained through the replica 
construction and the latter one occurs because the endpoints of the intervals are on the 
real axis. Indeed, the complex equations (73) and (74), which define the Riemann surfaces 
lZ n and 7 Z n , are invariant under complex conjugation. 

In [45, 46] the symplectic matrices which implement the dihedral symmetry of lZ n 
have been written explicitly and in App. C.3 this analysis has been extended to lZ n 
as well (the symmetry Z 2 is different in the two cases). These transformations act 
on the period matrix and reshuffle the characteristics, but the functions and f2 n [e] in 
(63) remain invariant. Moreover, both the transformations associated to the dihedral 
symmetry leave the coefficient r n (S) in (56) invariant. Thus, the terms in the sum (55) 
whose characteristics are related by one of these modular transformations are equal and 
the sum can be written in a simpler form by choosing a representative term for each 
equivalence class, whose coefficient is given by (56) multiplied by the number of terms of 
the equivalence class. 

Exploiting these symmetries, one can write the explicit expressions given in Sec. 6.3 
of [59] for 2 ^ n ^ 5. Beside the goal of having more compact analytic expressions, the 
dihedral symmetry is very helpful also from the numerical point of view because it allows 
to reduce the exponentially large (in n) number of terms in (55). 

Looking at Eq. (13) on the lattice, the Z n symmetry corresponds to the cyclic 
permutation of the n factors within each trace. Instead, the Z 2 symmetry comes from the 
fact that p 0 and p\ are not separately hermitian but the hermitian conjugation exchange 
them, so that pj is hermitian. However, as already noticed, such exchange leaves any 
term of the sum unchanged. 

3.3. Self-dual boson 

In this subsection we show that the expression (55) for Tr (pj)" of the Dirac free fermion 
is equal to the one for the compactified boson at its self-dual radius. 

The analytic formula for Tr (pj)” of the compactified boson for a generic value of 
the compactification radius has been derived in [ 6 ] by studying the partition function of 
the model on the Riemann surface lZ n . At the self-dual radius, it becomes (see Eq. (146) 


of [ 6 ] for 77 = 1 ) 



(57) 


where the matrices occurring in this expression have been defined in (50). The Riemann 
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0 


i x 77 
77 ix 



(2f) 


e 


£ 

0 


(2f)t 


(58) 


where in the first step we have used (3.6b) of [57] and in the second one —2r* = 2r — 2Q. 
Then, by specialising the addition formula reported in [55] (pag. 4) to our case, we find 


0 


i X 77 
77 iX 


1 

2 n ~ 1 


£(- 1) 4 * 4 


e 27rieQe 0 



2^lX>n(<S)0 

5 


0 

5 


(r) 2 , (59) 


where 


^W = ^(-l) 4 £ ' < 5 e 2 " i£ - s ' 
£ 



(60) 


In App. D we show that m n {8) can be written as the partition function of a classical Ising 
spin system, where £ play the role of the spin variables and the 8 are the local magnetic 
fields. In the same appendix we also employ standard transfer matrix techniques to prove 
that m n (8 ) = r n (8) (see (56) and (60)). 


3-4- Numerical checks 

In [59] we have already shown that Tr(p^ 2 ) n for n = 3,4,5 converges in the continuum 
limit to the CFT predictions (55). I 11 this paper we have given a set of more stringent 
relations (52) between each term in the sum for Tr(pJ) n appearing both in CFT and on 
the lattice. The goal of this subsection is to provide explicit numerical evidence of this 
term-by-term correspondence for n = 2,3,4. 

In order to evaluate numerically the traces of product of these matrices, we employ 
the techniques first developed in [60] for Trp^ and recently used to compute Tr(pJ) n in 
[59]. Indeed, being the tight-binding Hamiltonian (1) quadratic in the fermionic operators, 
the ground state reduced density matrix pa is Gaussian. Moreover, pA in (12) is Gaussian 
as well and, since the string operator Pc can be written as the exponential of a quadratic 
operator, the density matrix Pa 2 PaPa 2 is also Gaussian. Nevertheless, the sum of these 
two matrices in (12) is not Gaussian (this is indeed the main difficulty compared to bosonic 
models in which the partial transpose is itself Gaussian [61, 62, 63, 64]). By exploiting 
the fact that for Gaussian states all the information of the system is encoded in the 
correlation matrices, the computations can be performed in a polynomial time in terms 
of the total size of the subsystem. I 11 particular, in our case the correlation matrices of 
Pa and Pa 2 PaPa 2 can be obtained from the one of p^, as described in [26, 59]. 

The lattice computations have been performed in an infinite chain. The disjoint 
blocks A\ and A 2 have been taken with the same size i\ — £2 = 7) while the size of the 
block Bi separating them is d. Thus, the four point ratio (43) becomes 


x = 


e + d 


(61) 















J^/J 3 J 2 lat /J 2 
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Figure 3. The terms occurring in Tr (pf^) n for the free fermion (see (55)), according to 
the correspondence (52). Here we show n = 2 (top panels) and n = 3 (bottom panels). 
For each group of identical terms, only one representative has been plotted. In the left 
panels, the term with p = 0 has been divided by its CFT counterpart (5 = 0), in order 
to simplify the residual dependence on t\ and li- The extrapolated points (red crosses) 
are obtained through a fit of the data according to the scaling function (64) and they 
agree with the CFT predictions (solid lines). 


and configurations with the same value of l/d correspond to the same x. 
Referring to Eq. (14), let us introduce the following lattice quantities 


J? = Trp£, 


^[P? = 


Tr 


Tr 


71—1 


Po 


n Ppk 


k =1 


(62) 


which can be evaluated as explained in [59]. We also introduce their CFT continuum 
limit: 


J — c 2 

°n — c n 


X 


4t2 


2A„ 


O r-1 _ 0 [ e ](^)) 

■ = e(f W) ' (63) 

These CFT values are approached by taking configurations with increasing £, keeping the 
ratio £/d fixed. As discussed in Sec. 3.2.1, many terms in the sum (14) are equal because 
of the properties of the trace (in the continuum, this degeneracy is due to the dihedral 
symmetry of the Riemann surface). 

In order to deal with the finite size effects, we perform an accurate scaling analysis, 
as done in [60] for Trp^ and in [59] for Tr(pJ) n . From general CFT arguments it has been 
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Figure 4. The terms occurring in Tr(//J) 4 for the free fermion (see (55)), according 
to the correspondence (52). For each of the three groups of identical terms, only one 
representative has been plotted. In the top panel, the term with p = 0 has been divided 
by its CFT counterpart (5 = 0), in order to simplify the residual dependence on l\ and 
£ 2 - The extrapolated points (red crosses) are obtained through a fit of the data according 
to the scaling function (64) and they agree with the CFT predictions (solid lines). 


shown that these quantities display some unusual corrections to the scaling in £ described 
by a power law term with exponent 5 n = 2 h/n, being h the smallest scaling dimension of 
a relevant operator inserted at the branch points [65, 66, 67, 68]. For the Dirac fermion 
h = 1 and terms of the form £~ 2m / n are present, for any positive integer m. Because of 
the slow convergence of these terms (which becomes slower and slower for increasing n), 
typically it is necessary to include in the scaling function many of them. The most general 
f ini te-/ ansatz for f2 n takes the following form 



^[e] + 




p/ 


x) cai 2) 
— + 


£4/ 


\ (3) 

X) UJn 


'x) 


p/ r , 


+ . . . , 


(64) 


where fi[e] is defined in (63) and p and e are related through (53) and (54). For the 
J^ at /J n , a scaling function similar to (64) can be studied. Fitting the data with (64), 
the more terms we include, the more precise the fit could be. Nevertheless, since we 
have access to limited values of £, by using too many terms overfitting problems may be 
encountered, which lead to very unstable results. The number of terms to be included in 
(64) has been chosen in order to get stable fits. We find that every term D n [e] follows the 
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scaling (64) and the extrapolated value agrees with the corresponding CFT result. 

Our numerical results are shown in Fig. 3 for n — 2 (top panels) and n — 3 (bottom 
panels), while Fig. 4 is about the n — 4 case. As for the prefactor, the ratio J^ at / J n has 
been considered in order to eliminate the trivial dependence on t which survives in the 
continuum limit. The solid lines are the CFT predictions, which are given by (52). 

4. Conclusions 

In this manuscript we studied the moments of the partial transpose of the reduced density 
matrix Tr(pJ)" for two disjoint intervals in the conformal field theory of the massless Dirac 
fermion. Our main result is a closed analytic form for these moments of arbitrary order, 
i.e. Eq. (55). For n = 3,4, 5 this formula was anticipated in Ref. [59], but we extend here 
to arbitrary n and provide its full derivation. The analytic computation of the logarithmic 
negativity £ through the replica limit of (55) for even n e —> 1 is beyond our knowledge. 

It turned out that these moments are identical, for arbitrary order, to those of 
the compactified boson at the self-dual point. This equality comes from the explicit 
computation and we miss a proper understanding of this fact. It was already noticed that 
for the moments of the reduced density matrix of two disjoint intervals Trpjj, the result 
for the free fermion [41, 69] and the one for the compactified boson at the self-dual radius 
[38] are equal and very easy (they are both given by (42) with J r n (x) = 1). This is not the 
case for three or more disjoint intervals [45, 46]. This unexpected equivalence has been 
investigated in [45], where also other results have been found, based on the fact that r(x) 
is purely imaginary when x G (0,1). For the partial transpose, the period matrix t(x) of 
7 Z n in (50) has a non vanishing real part. Nevertheless, here we have shown that Tr ( p T ff )” 
for the free fermion is equal to the one for the self-dual boson, a property that does not 
follow from the analysis of [45] . The equality of all the moments obviously implies also the 
equality of the negativities. Since the negativity is directly measurable by means of tensor 
network algorithms (as e.g. done in [21, 7] for the Ising model), it would be very interesting 
to check numerically the identity between the negativity of the tight-binding model and 
the isotropic Heisenberg antiferromagnet (whose continuum limit is the self-dual boson). 
This is a highly non trivial prediction. 

We point out that an interesting technical byproduct of this paper is a one-to-one 
correspondence between each of the 2 n_1 terms appearing in the lattice formulation 
of the moments Tr (p 7 ^) n and the partition function of the free fermion on 1Z n with 
a particular assignment of fermionic boundary conditions. This correspondence has 
been explicitly checked against lattice numerical computations extending the numerical 
analysis of Ref. [59] where only the overall sum was considered. The consequence of 
this correspondence for the moments of both the reduce density matrix and its partial 
transpose in spin models with a free fermionic representation will be explored elsewhere 
[70]. 
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Appendices 

A. Reduced density matrix and its partial transpose on the lattice 

In this appendix we briefly review the main tools employed in this manuscript to study 
the partial transpose for free fermions on the lattice. 

Given the reduced density matrix (7), it is convenient to distinguish the terms having 
an even or odd number of fermionic operators in A 2 (notice that the parity of operator in 
A2 is the same of the operators in Af) by introducing 

Peven = ^1+^2 X^ Wl2 P° dd = X^ Wl2 O^ 2 ' ( 65 ) 

even odd 

Thus pa = Peven+Podd- The partial transposition with respect to A 2 in (10) acts differently 
on the two operators in (65). In particular [26] 

Peven = ^(-l)^ 2 ^ O x 0 2 , P^dd = 0 1 0 2- ( 66 ) 

even odd 

Defining the following Gaussian matrix pA (which is not a density matrix, being not 


hermitian) 

PA= 2 h+h O1O2, 

(67) 

and also p even 

and p 0 dd as done in (65) for p^, the partial transpose of Pa becomes 



pA Peven 1 Podd • 

(68) 


The matrices p eve n and Podd in ( 68 ) can be written through the string Pa 2 of the Majorana 
operator along A 2 . Since Pa 2 o x Q v Pa 2 = (—l )^' sA 2 a x -' y , one finds 

Peven = ^ (Aa T" PA 2 P aPA^) i Podd = ^ {ft A Pa 2 PaPa 2 ) • 

Thus, plugging (69) into (68), the final expression for (7} in (12) is obtained [26]. 


( 69 ) 
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In this appendix, by employing the formalism described in Sec. 2.2 and Sec. 2.3, we check 
the standard relation between the reduced density matrix of two disjoint intervals and its 
partial transpose 

Tr (pj)' 2 = Tr (p 2 „) , (70) 

From (12), it is immediate to observe that Tr (pj) 2 = Tr (paPa 2 PaPa 2 )- Then, by using 
(18), we can write 

Tr (paPa 2 PaPa 2 ) = (71) 

= j DCDCDp*Dpe^ cc e~P v (~Ci, ~C2 \pa\vu V2) (Vi, V2 \Pa 2 PaPa 2 \(i, C 2 ) 

= f DC DC Dp*Dpe~ cc e~ ,1 * v (~C U -C 2 IPAIP 1 , V 2 ) {vi, ~V2 \Pa\Ci, -C 2 ) 

= f DCD(D V *D V e-P<e-^(- G, V * 2 \p A \m, C 2 ) (Vu -C 2 MC 1 , V* 2 ) 

= I DCDCD7fD V e-^ + ^e-^ + ^ <-Ci, mMvi, C 2 ) (hi, -C 2 MC 1 , V 2 ), 

where in the last step the change of variables C 2 ~t C 2 , h 2 ~t pi, has been employed. Then, 
by noticing that the relations in (18) can be slightly modified as follows 

-1= |dC*dCe rc |C){-CI. — Tr O — |d<*d<e« ((| 6 |(), (72) 

we can conclude that (72) is exactly Tr (p\)- 

C. On the Riemann surface lZ n 

In this appendix we derive with all the details some results about the Riemann surface 
lZ n and its period matrix r(x) in (50) that are employed in the main text. 

Given the two disjoint intervals A\ = {u\,V\) and A 2 = ( u 2 ,v 2 ) whose endpoints are 
ordered as U\ < V\ < u 2 < v 2 , Tr p n A is the partition function of the model on the Riemann 
surface 7 Z n which is dehned by the following algebraic curve in C 2 (parameterised by the 
complex variables z and y) [38] 

y n = (z - ui)(z - u 2 ) [{z - Vi)(z - v 2 )] n l . (73) 

The Riemann surface lZ n , which is an n sheeted cover of the complex plane, has genus 
7 i—l and it has been studied in detail in [71, 46], where its generalisation to any number 
of disjoint intervals (whose genus is (n — 1)(N — 1) for N intervals) has been also discussed. 
By going around u\ and u 2 clockwise, one goes from the j-th to the (j + l)-th sheet, while 
going around V\ and v 2 clockwise, one moves to the (j — l)-th one. 

In order to compute Tr(pj) n , one has to find the partition function of the model on 
a different Riemann surface 7 Z n , which is obtained by exchanging u 2 <H- v 2 in (73), i.e. 

y n = (z - u 1 )(z - v 2 )[[z - Vi)(z - u 2 )] n 1 


(74) 
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The Riemann surface lZ n has genus n — 1 and the sheets are joined in a different way with 
respect to 7Z n . Indeed, by encircling w 2 clockwise we move from the j'-th to the (j — l)-th 
sheet while by encircling u 2 clockwise the (j + l)-th sheet is reached. 

While the period matrix t(x) of 7 Z n is purely imaginary (see (45)), the period matrix 
f{x) for 7 Zn has a non vanishing real part. In Sec. C.l we show that Re[f(a;)] has a 
very simple form. In Sec. C.2 we consider the moments Tr (pj)" in the regime of small 
intervals x —> 0 and in Sec. C.3 we provide a detailed discussion of the Z 2 part of the 
dihedral symmetry for 7 Z n . 


C.l. The real and imaginary part of the period matrix r(x) 


In this subsection we want to write explicitly the real and the imaginary part of the period 
matrix f(x) given by (50). The real part Re[r(a;)] turns out to be a simple tridiagonal 
matrix with half-integer entries. 

Let us introduce the following ratios of hypergeometric functions, which enter in the 
expressions for the period matrices r{x) and f(x) (see (45) and (50)) 


/ \ _■ 2 -Pi(r, 1 -r; T,l~x) 

tax) = i 


T r (x) = T r 


X 


2 Fi(r, 1 -r] l;x) ’ '\x-l 

where 0 < r < 1 and x G (0,1). Moreover, we also define 
= T(l-2r) 2 Fi(r,r;2r; 1 - x) 


= a r (x) + i j3 r (x), (75) 


r(i 


2 Fi(r,r; l;x) 


« M = r(2r-l) _ v _ 2r 2 Fi(l-r,l-r;2(l-r);l-g) 

rU _ T(r) 2 1 j 2 Fi(r, r; 1; x) 

By employing the expressions given in (87) of [6], one finds that 


(76) 

(77) 


a r (x) = sin(7rr) [A r (x) + B r (x )\, /3 r (x) = cos(nr)[A r (x) — B r (x)]. (78) 


At this point we need the following identity (see e.g. Eq. (1) at pag. 108 of Ref. [72]) 
2 Fi(a, b ; c; z) = ^ ^ —jl 2-^1 (a, 6 ; a + b - c + 1 ; 1 - z) 


T(c — a) T(c — b) 
T(c) T(a + 6 — c) 


(79) 


(1 - z) 


c—a—b 


<Fi(c — a, c — 'b; c — a — b T 1; 1 — z), 


T(a) T(b) 

which holds for |arg(l — z)\ <7 r. By specialising (79) to the case of (a, b, c ) = (r, r, 1) and 
z = x G (0,1), from (76) and (77) one finds that 


A r {x) + B r {x) = 1. 


(80) 


From this expression it is clear that the x dependence disappears from the real part of 
f r (x) and hence from the period matrix. Indeed, using (78) and (80) in (45) one gets (see 
also (50)) 

2 1 

f(x) itj = -^2sm(irk/n)? k/n (x) cos[2n(k/n)(i - j)} = - [Q].,. + i [l(x)\.., (81) 

1 k =1 
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where the sum giving the real part can be explicitly performed, finding that the matrix 
Q has integer elements which read 


[2]*J — 2Si,j j|,l> 


(82) 


namely Q is a symmetric tridiagonal matrix, with 2 on the main diagonal and —1 on the 
first diagonals. On the other hand, the imaginary part can be written as follows, 


n— 1 


[lW]j, y = - ^sin(rt/n) A/„(i)cos[2jr(fc/n)(; - j)]. 


(83) 


fc=l 


with 


Pr(x) = cos(vrr), f r (x) = ^^ x r 2 F 1 (r, r; 2r; x). (84) 


(1 - x) r 2 Fi(r, r; T,x) 


r(l - r) 


This expression for T{x) can be further simplihed. Plugging (80) into the second 
expression of (78), we have 


/ 3 r (x ) = cos(7rr) [2 A r (x) — l]. 

For 0 < x < 1 we can rewrite A r (x) as follows [72] 

1 


A r (x) = 


-7-\" 

2 cos(7rr) 


—inr I inr -r 2^1 Ob V, 1] l/x) 

c- | o iL> 


2 Fi(r,r\ l;i) 


(85) 


( 86 ) 


Since A r (x) is real for 0 < x < 1, the vanishing of its imaginary part gives 


x' 


tan(7rr) Re 


2 Fi(r,r; 1; l/x) 
2 Fi(r,r; l;x) 


+ Im 


2 Fi(r,r; 1; l/x) 
2 Fi(r,r; T,x) 


= tan(7rr). (87) 


On the other hand, by writing A r (x) as its real part and plugging the resulting expression 
in (85), one finds 


A-0*0 = — 

x r 


cos(7rr) Re 


2 Fi(r,r; 1; l/x) 


2 F 1 (r,r; l;x) 
Finally, using (87) we can write 


— sin(7rr) Im 


2 Fi(r,r; 1; l/x) 
2 Fi (r,r; l\x) 


( 88 ) 


f3 r (x) = cos(7rr) 


2 Fi(r,r; 1; l/x) 
sin(7rr) 0 2 Xi(r, r; 1; x) 


Im 


(89) 


The matrix Z(x) can be easily written by plugging (89) into (83) and noticing that the 
sum over the cosine vanishes. The result reads 

[l(x)\. . = - - V x- k/n cos[2t x{k/n){i-j)\. (90) 

L h,j n l p (k/n, k/n: 1; x) J 


k =1 


The result (82) is employed in Sec. 3.3 and in App. C.3. 
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C.2. Short intervals regime 


In this appendix we study the Tr ( p 1 /)" for the free fermion (55) in the limit of short 
intervals, i.e. when x —> 0. 

In the expression (55), only Riemann theta functions with e = 0 occur, which are 
given by 


0 


(r(a;)) =1+^2 

m/O 


gi7r(m- Q- m-\-26 m) n m- Z- m 


(91) 


where Q is independent of x. Expanding 6k/n( x ) in (84) for x —* 0, one finds 


Pq(x) = _. 8m ^ rg ) [log(x) + 2y E + 'f’iq) +^(1 “?)] +0(x). (92) 


Plugging this expansion into (83) and (91), one gets that the leading term is x rn Q rn . 
The exponent m ■ Q ■ m for m G Z n_1 has been already analyzed in [39], finding that its 
minimum is 1, which is obtained for the following vectors 


_ h _ 

m± = ( 0,..., 0,±1,...,±1, 0,..., 0), m± ^ 0, 0 ^ j x < j 2 ^ n - 1, (93) 

h 

namely m± ■ Q ■ m± = 1. Then, by applying again the results of [39] (notice that the 
vectors m + and m give the same contribution) to (91), we find 


0 




■i) 


2 <5m+ 


m+ 


sm (n(j 2 - ji)/n) 


+ • • •, 


(94) 


where the dots denote o(x) terms. Thus, for the generic term occurring in the sum (55) 
we have 



x ^ 1 + (-1) 2<5 '™+ 

2 ” 2 Sil ' 2 ( 7r (i 2 - jl)/n) + " ' ’ 



(95) 


Plugging this result into (55), we get the first term of the x —» 0 expansion of Tr (p^ 2 )”, 

which is 


Tr (pj)” = j n 


1 - 


X 


2”-’ n 2 


Y, 


i + (-i) 


_i \2S m+ 




sin (7r(j2 ~ ji)/n) 


+ 


(96) 


C.3. The Z 2 part of the dihedral symmetry of lZ n 

In this subsection we briefly discuss the most peculiar aspect of the dihedral symmetry 
for the Riemann surface 7 Z n occurring in the computation of Tr (p^ 2 )™ (see Sec. 3.2.1). 

The Riemann surface 7 Z n has a dihedral symmetry Z n x Z 2 due to the invariance 
under cyclic permutation of the sheets (Z n ) and the complex conjugation (Z 2 ). For 
a genus g Riemann surface, the modular transformations are given by the symplectic 
matrices Sp(2g,Z) [56, 57]. The dihedral symmetry can be identified with a subgroup 
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of the modular transformations acting on 7 Z n which has been discussed in [45, 46]. In 
particular, these peculiar modular transformations map the a cycles among themselves 
and the b cycles among themselves, leaving the period matrix t(x) unchanged. 

Also the surface lZ n has a dihedral symmetry Z n xZ 2 but, while the cyclic permutation 
(Z n ) is exactly the same one discussed above for lZ n , the complex conjugation is slightly 
different because it mixes a and b cycles. Let us remind that the complex conjugation 
corresponds to reverse the order of the sheets and to reflect all of them with respect to 
the real axis. Considering the canonical homology basis {a r , 6 r ; 1 ^ r ^ n — 1} introduced 
in Sec. 3.2 (see Fig. 2) [46], this transformation acts as follows 

= M inv • ( ~ ) , M inv = ^ g 1 | 0 j e Sp(2(n - 1), Z), (97) 



where we introduced the notation of the double-headed arrow above a matrix to indicate 
that the columns have to be taken in the reversed order (I*, is the k x k identity matrix 
and Q is given by (82)). Under the symplectic transformation defined in (97), the period 
matrix f (50) changes as follows 

t' — I • t • I — Q • 1 — t — Q — —77 + iZ = — t*, (98) 


while for the characteristics of the Ricmann theta functions we have 


Notice that the powers of M inv read 


I 0 

•fA -fA 

Q I 


M 2k ~ l = 


inv 


•fA 

I 0 

•h- e I , 


(99) 


M 2k — I ^ ^ 

inv 1 -2 kQ I I ’ 


-(2k-l)Q I 

(in particular, notice that M 2 nv ^ I) so that, by applying (98) k times one finds 


( 100 ) 


f (fc) = t — k Q. 


( 101 ) 


As for the change of the Riemann theta function under the modular transformation in 
(97), because of the particular form of Q, it is easy to show that for k even it is left 
invariant (a part for an overall sign), while for k odd it becomes its complex conjugate, 
up to an overall sign. Since in our formulas the modulus of the Riemann theta function 
always occurs, the terms occurring in our sum over the characteristics are invariant under 
this transformation. Thus, Mj nv can be the modular transformation representing the Z 2 
of the dihedral symmetry, even if Mf nv ^ I. 
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In this appendix we show the equality m n (S) = r n (d) between the coefficient m n (S) in 
(60), coming from the self-dual boson approach, and the coefficient r n (S) in (56) occurring 
in the expression obtained through the free fermion analysis. 

Considering m n (S) in (60), the expression in the exponent can be written as follows 


^ n —2 n —1 n—1 

£■ — ■£-\-£-S = — 'y ' £i £j+1 + y ' £} + y ' £iSi. (102) 

2=1 2=1 2=1 

Then, defining the spin variables <7j = 4£j — 1 = ±1 and local magnetic helds hi = 4 Si, 
we hnd that m n (S ) is equal to the partition function Z n of n — 1 Ising spins in a binary 
magnetic held hi = 0, 2 (a part for the first and last site), which reads 


= e T(E, 


i+n) 


exp 


7r 

l — 

4 


n—2 n—1 

— CTjCTj+i + CTihi + cq + cr n _i 


2=1 


2=1 


(103) 


Given this Ising model representation, we can compute m n (S) through standard transfer 
matrix techniques. Following [73] , let us introduce the conditioned partition function with 
the last spin set to /i, namely 


Z n (n) = e 4 &i hi+n ) y>xp 


17T 

J 


(105) 


n—3 n—2 

■ ^ ^ <7?: CTj_|_1 — flG n -2 + ^ CTj/ij + /l h n _i + (Ti + /O 
2=1 2=1 

(104) 

Then, by adding one spin e to the partition function, it becomes 
Z n+1 (e)=e?( E *H x 

n—3 n—2 

— ^ ~] (TjCTj-i-i — /aa n -2 — /i e + ^ <J % h{ + // /r n _i + e /i re + ay + e 
2=1 2=1 

After some algebra, one realizes that 

Zn+i(+) = (~l) 2Sn [Z n (+) - Z n (-)} , z n+1 (~) = Z n (+) + Z n (~). 

and these relations give the transfer matrix 


x 5Z ex p 

H=± 1 <T 


17T 

T 


(106) 


= * ( f - 1 ) 25 - f - 1 ) 25 

1 ^ y/2 \ 1 1 


<5 6 { 0 , 1 / 2 }. 


We also need the conditioned partition functions for a single spin, which read 

^(+) = - e 2 hl , Z 2 (—) = 1. 

Then, the partition function for n — 1 spins in (103) is given by 
Z n = Z n (+) + Z n {—) = 

" ( - 1 1) ”‘) = ^ (1 1) ri 

k=n —1 x ' k=n —1 


(107) 

(108) 
(109) 


a?- 1 (11) n Th) 
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In order to compute the matrix product in (109), it is convenient to perform a change of 
basis which diagonalises T(0), namely 


T(5) = U ] T(S) U, 


U 


1 /l + i 

2 1 1 - i 


-(l + i) 
1 -i 



( 110 ) 


From (107) and (110) we can explicitly write the transfer matrix in the new basis 


T( 0 ) 


J_ /l + i 0 \ 

\/2 ( 0 1-ij 


n 1/2) 


±( °. 
\/2 yl + i 



and therefore the partition function (109) can be rewritten as follows 

i 

Z n = 2?- 1 (l 1)T(0) J] T(S k ) 

k=n —1 



an) 


( 112 ) 


Now, it is convenient to move all T(0)’s in the product of (112) to the left of all the 
T(l/2)’s. To do this, one observes that T(0)T(l/2) = T(1/2)T(0) _1 and T(1/2)T(0) = 
T(0) _1 T(l/2). By employing the latter rule, one finds 

Z n = 2?- 1 ( 1 1 ) T(0) 1+E > (1 “ 25i)( " 1)Ej=i T(l/2) E ‘ 2Si f | V (113) 


The factor 1 — 2 Si within the sum in the exponent of T(0) selects only the sites where 
Si = 0, while the other factor (—1)^'=* 2Sj counts all the T(l/2)’s on the left of site i. The 
exponent of T(0) can be rewritten as 


n— 1 

s{6) = 1 + ^(1 - 28i)(—l)^=i 26j 

i =1 


1 - (-l)Si25i 

1 + ----- 

2 


n— 1 

+ £(-i) ES ' 2Si - 

i= 1 


(114) 


Since T(0) is diagonal, its powers can be easily performed. Moreover, since T(l/2) 2 = I 2 , 
every integer power of T(l/2) is simply I 2 if the exponent is even, and T(l/2) otherwise. 
Thus, we have 


T(0) s(<5) 



(115) 


T( l/2) 2 ^ Si 


[1 + (_ 1 )Ei^l /2 e~'l[l- (—l)^ 2(S *]/2\ 

e T[l-(-!)£. 2^/2 [i + (_i)E t 25 s]/2 ) 


(116) 


Finally, (113) tells us that we just need to multiply this two matrices and to sum the four 
elements of the resulting matrix. After some of algebra, we get 


m n (S) = T' 2 


1 + ( —l)£i 2 <5,: 


COS | ^ s (<5) 


+ 


1 (l)£i 2 * 


cos 1 ^ (s(<5) 


1) 


• (117) 


By inspection of the two cases of JT 25i even or odd, it is clear that (117) equals (56). 
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